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I.  Introduction 


During  September  of  1980  and  January  of  1981  the  Naval 
Postgraduate  School  (NPS)  and  Aerovironment ,  Inc.  (AI)  conducted 
an  intensive  series  of  tracer  measurements  off  Ventura,  California 
in  the  Santa  Barbara  Channel.  The  main  purpose  of  these 
experiments  was  to  parameterize  Gaussian  and  trajectory  dispersion 
models  for  the  overwater  regime.  A  secondary  purpose  was  to  build 
a  data  base  for  future  model  development. 

The  experimental  procedure  was  as  follows:  SF^  gas  was 
released  from  the  NPS  research  ship  RV/Acania  at  a  rate  of 
approximately  50  lbs/hr.  The  ship  was  anchored  at  a  point 
approximately  5  miles  from  shore.  The  position  and  dimensions  of 
the  plume  were  determined  at  three  downwind  distances  by  analyzers 
aboard  an  aircraft  and  van,  and  by  collecting  gas  samples  at  fixed 
shore  locations  and  from  a  small  boat.  The  plume  analysis  and 
shore  based  meteorological  measurements  were  performed  by  AI .  A 
complete  characterization  of  the  overwater  meteorology,  including 
stability,  surface  fluxes,  and  horizontal  wind  direction  variance 
was  obtained  by  NPS . 

Data  volumes  for  both  experiments ( 1 )  and  a  final  report 
on  the  results  of  the  analyses  will  be  available.  The  purpose  of 
this  report  is  twofold:  to  present  the  wind  direction  variance 
results,  which  are  not  included  in  the  former  data  reports,  and  to 
present  results  comparing  measures  of  overwater  stability  and  the 
wind  variance.  In  order  to  describe  the  conditions  for  which  the 
data  were  obtained,  a  brief  section  describing  the  locality  and 
observed  meteorological  conditions  is  included. 
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II.  Meteorological  Conditions  during  the  Experiment 

The  area  near  Ventura,  where  the  tracer  experiments  were 
performed,  is  representative  of  the  general  channel  islands  region 
but  not  of  the  majority  of  the  California  coastal  region.  It  is 
an  area  where  the  weather  is  controlled  by  the  interplay  of 
numerous  local  influences*  and  the  mesoscale  features  which  are 
typical  of  the  rest  of  the  coast.  The  area  lies  within  the 
embayment  formed  by  the  Santa  Barbara  Channel,  Santa  Monica  Bay, 
and  the  Gulf  of  Santa  Catalina.  The  Oxnard  plane  is  approximately 
20  miles  long,  extends  from  the  beach  to  five  to  ten  miles  inland 
and  is  surrounded  by  coastal  hills  with  peaks  two  to  three 
thousand  feet  high.  Within  40  miles  inland  mountains  extend  to 
8, 000  feet . 

The  geographic  features  cause  many  effects,  the  major  ones 
being: 

1.  The  mountains  to  the  North  act  as  a  partial  barrier  to 
air  mass  movement  from  the  North. 

2.  The  mountains  and  the  general  east-west  orientation  of 
the  shore  produce  a  complex  pattern  of  eddies. 

3.  Surrounding  high  hills  contain  the  cool  moist  marine  air. 
Only  infrequent,  strong,  synoptic  air  mass  changes 
displace  the  marine  air. 

4.  The  surrounding  hills  and  mountain  valleys  channel  the 
air  flow  into  complex  patterns. 

*  Much  of  the  description  of  the  local  meteorology  presented  here 
is  exerpted  from  Pacific  Missile  Test  Center  documents,  which  are 


listed  in  the  reference  section  as  reference  2. 
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Heating  of  local  slopes  "draws  in"  the  afternoon  sea 
breeze,  and  nightime  cooling  and  the  resulting  drainage 
reverses  the  local  flow.  The  diurnal  cycle  is  very 
persistant . 

6.  During  periods  of  storms  with  strong  southerlies,  the 

common  rain  pattern,  the  wind  is  turned  to  southeasterly 
and  accelerated. 

Because  of  the  complexity  of  the  local  influences  any 
measurements  made  in  this  area  are  not  representative  of  the  open 
coast  but  are  location  specific. 

It  is  difficult  to  obtain  sufficient  data  in  a  given  area  to 
parameterize  a  model  for  all  conditions  that  may  be  encountered. 
This  experiment  was  conducted  for  two  approximately  two-week 
periods,  one  in  the  late  summer  and  one  in  the  winter.  With  a 
limited  data  set  one  must  compare  the  meteorological  conditions 
encountered  with  typical  conditions  for  the  locality  then  attempt 
to  extend  the  results  to  other  conditions  that  are  frequently 
encountered.  Extreme  conditions  cannot  normally  be  treated.  For 
comparison  purposes,  we  list  below  the  normal  seasonal  conditions 
in  the  Ventura  area  and  conditions  encountered  during  the 
experiments.  This  will  constitute  the  meteorological  framework 
for  discussing  the  tracer  experimental  results. 
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A .  Seasonal  Conditions  in  the  Ventura  Area 
Summer 

The  North  Pacific  Semipermanent  Subtropical  High  lies  to  the 
west  of  the  area  and  controls  the  synoptic  scale.  Clockwise  flow 
around  the  high  produces  northwester lies  along  much  of  the  coast; 
this  flow  is  turned  more  westerly  by  the  Point  Conception/mountain 
geographic  influence.  The  local  sea  breeze  turns  the  wind  further 
to  a  westerly  or  southwesterly.  The  general  onshore  flow  is  aided 
by  the  inland  thermal  trough  which  is  created  by  overland  heating. 
Strong  subsidence  creates  the  prevalent  capping  inversion.  The 
occasional  passage  of  weak  upper  level  troughes  will  dissipate  or 
lift  the  inversion  for  short  periods  of  time. 

Fall 

The  building  of  high  pressure  in  the  Great  Basin  causes 
frequent  Santa  Ana  conditions.  The  pattern  of  storms  and  upper 
level  westerlies  moves  further  south  breaking  up  the  summer 
pattern.  Frontal  passage  becomes  more  frequent  and  the 
subtropical  high  becomes  displaced  or  shrinks,  resulting  in  a 
break  up  of  the  marine  inversion. 

Winter 

Frontal  passage  becomes  much  more  frequent  and  strong  surface 
westerlies  often  follow  the  passage.  Santa  Ana  winds  can  still 
occur  when  the  surface  pressure  in  the  Great  Basin  becomes 
sufficiently  high.  Also,  the  Pacific  High  and  capping  inversion 
can  reform  between  frontal  passage  occurrences. 
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Spring 


As  the  storm  pattern  moves  north,  the  Pacific  High  again 
becomes  the  dominant  feature.  Cold  lows  pass  frequently,  followed 
by  strong  westerlies. 

B .  Wind  Climatology 

Although  the  strong  local  influences  make  it  difficult  to 
accurately  predict  the  wind  in  the  Ventura  area,  the  seasonal 
patterns  are  quite  r eproducable .  This  makes  a  climatology  useful 
for  assessing  expected  onshore  impact.  Only  a  wind  climatology  is 
presented  here  since  that  is  our  primary  concern. 

Surface  wind  roses  for  each  month  are  shown  in  Figure  1. 

Wind  speed  is  indicated  by  the  width  of  each  "vector",  wind 
direction  by  the  angle,  and  frequency  of  occurrence  by  the  length. 
The  numbers  in  each  wind  rose  circle  are  the  percent  of  the  time 
the  wind  is  <4  knots. 

Average  prevalent  wind  as  a  function  of  time  of  day  and  time  of 
year  are  shown  in  Figure  2.  The  dashed  lines  show  wind  direction, 
where  a  vector  from  top  to  bottom  on  the  page  indicates  a  north 
wind.  The  solid  lines  are  wind  speed  isopleths.  The  heavy  dashed 
lines  approximately  divide  the  land  and  sea  breeze  regimes. 

Of  special  interest  are  the  frequency  of  strong  winds,  the 
frequency  of  occurrence  of  Santa  Anas,  and  the  most  frequently 
observed  direction.  This  information  is  given  in  Table  1. 
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C .  Weather  Description  for  the  Test  Periods 

The  following  information  is  derived  from  the  daily  weather 
maps,  weekly  series,  published  by  the  National  Oceanic  and 
Atmospheric  Administration.  The  periods  of  interest  are 
9/24-9/29,  1980  and  1/6-1/15,  1981.  For  completeness  we  describe 
the  situations  for  periods  slightly  longer  than  the  duration  of 
the  experiments. 

9/24-9/29:  Throughout  the  period  a  strong  thermal  low  persisted 

in  northwestern  Mexico. 

9/22:  Surface  highs  were  centered  over  the  Northwest  and 
Wyoming  and  pressure  gradients  in  the  western  US  were  weak. 

Strong  onshore  flow  was  not  expected. 

9/23:  Surface  highs  were  located  over  Colorado.  Pressure 
gradients  were  even  weaker. 

9/24:  Surface  highs  were  moving  in  off  Canadian  coast.  A 
weak  upper  level  ridge  was  forming  in  the  same  area  and  a  weak 
trough  forming  over  north  central  US.  Surface  gradients  were 
increasing  off  the  California  coast  causing  weak  onshore  flow. 

9/25:  Surface  highs  were  located  over  the  northern  portion  of 
the  continent.  The  upper  level  ridge  was  strengthening  but  only 
over  the  area  north  of  California.  The  trough  over  central  US  was 
strengthening.  Only  very  weak  surface  gradients  existed  in  the 
coastal  region. 

9/26:  Surface  low  was  on  the  Washington  coast  and  very  large 


high  covering  the  Midwest.  The  upper  level  ridge  over  western 
Canada  was  strong  but  did  not  extend  south  to  California;  the 
trough  was  moving  to  the  East  coast.  The  gradients  were  somewhat 


higher,  weak  onshore  flow  was  expected. 

9/27:  The  situation  was  much  the  same  as  on  9/26  with 
additional  highs  over  the  northern  US  and  a  weak  front  moving  in 
off  Washington. 

9/28:  The  pattern  continued  with  the  addition  of  a  weak 
frontal  passage  north  of  San  Francisco  and  some  trailing 
precipitation. 

9/29:  A  large  surface  high  was  centered  over  Idaho  and 
Nevada.  The  upper  level  North-South  ridge  moved  east  to  extend 
from  Canada  into  Montana.  Location  of  the  high  caused  weak  Santa 
Ana  conditions. 

9/30:  The  thermal  low  over  nothwestern  Mexico  weakened, 
essentially  disappearing.  The  whole  of  western  US  was  dominated 
by  surface  high  pressure. 

General  Comments:  The  whole  period  was  dominated  by  the 
Pacific  high.  All  frontal  activity  was  to  the  north  of 
California.  The  thermal  low  over  Mexico  was  not  strong  enough  to 
produce  a  dominant  onshore  flow,  the  surface  pressure  gradients  in 
the  coastal  region  being  very  weak.  The  local  flow  was  dominated 
by  the  diurnal  land-sea  breeze  cycle.  Low  subsidence  inversions 
were  present  under  the  dominant  high  pressure.  Weak  Santa  Ana 
conditions  can  occur  under  these  conditions. 

1/6-1/15:  We  do  not  give  a  day  by  day  description  as  for  the 

previous  period.  The  two  periods  were  similar  and  general 
comments  will  suffice  to  describe  the  situation.  The  Pacific  high 
was  unusually  strong  producing  a  mini-drought  for  what  is  normally 
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the  beginning  of  the  rainy  season  for  California.  Frontal 
passages  were  again  far  to  the  north.  There  was  no  well 
established  onshore  flow  regime  except  for  a  short  period  during 
1/12-1/13  when  the  surface  gradient  in  the  Southern  California 
area  increased.  As  before,  the  land-sea  breeze  cycle  would 
dominate,  however,  fairly  persistent  highs  over  the  inland  western 


US  strengthened  the  offshore  flow  so  that  periods  of  sea  breeze 
were  shortened. 

Most  Frequent  Santa  Ana  Occurrance 


Month 

Direction 

Wind 

Speed 

% 

%  greater 
than  21  knots 

Average 

No .  days 

Maximum 
No .  days 

JAN 

NE 

10 

15.5 

2.3 

9.3 

16 

FEB 

W 

9 

12.3 

1.8 

5.2 

12 

MAR 

W 

10 

18.3 

1.1 

2.8 

8 

APR 

w 

10 

26.  7 

1.4 

0.6 

2 

MAY 

w 

9 

28.7 

0.3 

0.3 

2 

JUN 

w 

8 

27.5 

0.0 

0 

1 

JUL 

w 

7 

25.2 

0.0 

0 

0 

AUG 

w 

8 

23.6 

0.0 

0 

0 

SEP 

w 

7 

19.5 

0.1 

0.4 

4 

OCT 

w 

7 

16.3 

0.5 

2.7 

8 

NOV 

NE 

7 

12.1 

1.0 

7.0 

19 

DEC 

N 

5 

12.8 

1.4 

9.3 

18 

Table  1.  Monthly  averages  of  the  most  frequently  observed  wind  direction, 
percent  of  time  the  wind  speed  is  greater  than  21  knots,  and 
the  number  of  days  of  Santa  Ana  winds  per  month.  The 
maximum  number  is  the  maximum  observed  over  a  ten  year 
period. 
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III. 


Wind  Direction  Variance 


The  wind  direction  was  measured  by  a  vane  mounted  on  the 
RV/Acania  20m  above  the  water  line.  The  vane  responded  to  both 
the  wind  and  to  ship  roll?  since  the  ship  normally  pointed  into 
the  wind  during  the  measurements,  pitch  had  little  effect.  Ship 
direction,  measured  by  a  gyrocompass,  and  ship  roll  were 
determined  and  logged  by  a  computer  approximately  every  1.5  sec. 

A  pendulum  mounted  on  the  roll  axis  of  the  ship  was  used  to 
measure  roll.  The  roll  induced  transverse  component  of  wind  at 
the  vane  is  the  roll  rate  times  the  distance  from  the  roll  axis  to 
the  vane. 

The  wind  data  were  processed  and  stored  in  two  ways .  Mean 
wind  direction  and  speed  were  determined  for  half-hour  periods.  A 
correction  for  the  ship  direction  was  applied  to  the  wind  data 
when  determining  this  average.  A  roll  correction  was  not  made 
since  the  short  period  roll,  approximately  10-sec,  averages  to 
zero  over  a  half-hour  period.  Wind  direction  was  also  averaged 
for  10-sec  periods  and  these  averages  stored  for  later  processing 
to  determine  the  variance.  The  wind  data  was  corrected  for  the 
ship  roll  component  and  ship  heading  for  each  1.5-sec  wind 
determination . 

Hburly  averages  of  wind  direction  and  standard  deviations  are 
presented  in  Table  1.  Note  that  the  averages  are  not  exactly  one 
hour  since  they  are  compiled  from  the  half-hour  data  acquisition 
periods  and  there  is  dead  time  for  processing  and  printing  at  the 
end  of  each  acquisition  period. 


Table  2 


Hourly  averages  of  wind  directions 
and  their  standard  deviations. 


Average  Average 


Data /Time 

Direction 

Sigma 

Data/Time 

Direction 

Sigma 

9/24 

1/9 

1109-1205 

231 

13.3 

1119-1221 

278 

10.6 

1206-1301 

231 

12.9 

1239-1339 

273 

6.7 

1303-1357 

252 

11.5 

1339-1439 

281 

3.4 

1358-1453 

262 

7.7 

1440-1539 

279 

4.8 

1455-1549 

254 

8.0 

1540-1639 

276 

9.2 

1551-1645 

262 

7.1 

1640-1739 

272 

3.1 

1743-1837 

275 

6.0 

1/13 

9/27 

0822-0948 

11 

21.1 

0647-0743 

43 

6.3 

1019-1119 

299 

11.2 

0917-1013 

279 

10.2 

1209-1309 

298 

5.6 

1110-1206 

270 

6.6 

1310-1409 

287 

5.9 

1207-1302 

274 

5.6 

1409-1509 

271 

11.8 

1303-1358 

273 

4.7 

1509-1559 

255 

7.4 

1359-1454 

272 

4.0 

1559-1659 

240 

8.6 

1456-1550 

275 

2.4 

1551-1646 

273 

2.9 

1/14 

1647-1742 

268 

2.8 

1100-1200 

181 

46.5 

1774-1838 

269 

3.6 

1200-1300 

213 

23.6 

1300-1400 

194 

15.6 

9/28 

1400-1500 

246 

22.6 

1239-1335 

231 

10.3 

1340-1436 

247 

7.8 

1/15 

1437-1532 

257 

7.4 

1411-1500 

204 

17.3 

1534-1628 

253 

4.9 

1522-1622 

282 

32.7 

1712-1808 

260 

4.4 

1622-1722 

289 

6.5 

1809-1904 

257 

4.5 

1/16 

9/29 

0908-1008 

344 

12.7 

1149-1245 

229 

5.1 

1020-1120 

324 

5.6 

1251-1347 

253 

5.0 

1120-1220 

336 

6.2 

1349-1443 

262 

3.3 

1220-1320 

335 

8.9 

1444-1540 

261 

3.9 

1320-1420 

319 

17.7 

1558-1654 

266 

2.4 

1420-1520 

298 

4.6 

1656-1750 

268 

5.2 

1520-1620 

295 

4.3 

1751-1846 

285 

3.8 

1621-1720 

295 

4.8 

1721-1820 

294 

8.0 

1/6 

1821-1920 

299 

6.8 

1125-1225 

289 

13.7 

1921-2020 

271 

53.2 

1325-1425 

297 

11.6 

2021-2120 

195 

44.3 

1425-1542 

292 

22.6 

2121-2220 

277 

8.8 

1542-1642 

301 

13.3 

1642-1742 

279 

9.5 

1742-1842 

281 

17.1 

1842-1942 

342 

50.7 

12 


IV.  Stability  Classification 

It  is  common  to  use  some  measure  of  the  stability  of  the 
atmosphere  to  determine  the  expected  dispersion  of  a  plume.  The 
relationship  between  the  three  dimensional  spreading  and  stability 
is  complex  and  will  not  be  discussed  here.  We  present  the  various 
classification  schemes  and  their  interrelations.  Since  the 
Pasquill  classification  is  so  widely  used  we  will  show  the  method 
we  use  to  modify  it  for  the  overwater  regime.  Our  normal  method 
for  determining  stability  is  to  calculate  the  Monin-Obukhov 
length.  In  what  follows  we  will  compare  the  measured  wind 
direction  variance  to  this  length  and  to  the  modified  Pasquill 
classes  to  assess  the  usefulness  of  the  classification  schemes. 

Atmospheric  stability  depends  on  the  heat  flux,  which  is  a 
result  of  the  temperature  gradient,  and  the  amount  of  mechanical 
turbulence  which  is  generated  by  the  wind.  Several  schemes  for 
specifying  the  stability  are  in  current  user  all  require  a 
measurement  or  an  estimate  of  the  flux  and  turbulence.  The 
methods  we  briefly  describe  here  are  the  Monin-Obukhov  length,  the 
Richardson  number,  the  Pasquill  classifications,  and  an  EPA  method 
for  relating  the  Pasquill  class  to  the  temperature  gradient. 

A.  Monin-Obukhov  Length 

The  Monin-Obukhov  length  can  be  determined  from  measured 
gradients  of  wind,  temperature,  and  water  vapor  content  in  the 
surface  layer.  Gradients  are  difficult  to  measure,  especially 
from  a  moving  platform  such  as  a  ship.  Under  such  conditions 
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it  is  preferable  to  use  the  bulk  aerodynamic  method,  the 
measurements  required  being  sea-surface  temperature  and  wind 
speed,  temper ature ,  and  water  vapor  content  at  a  reference  height 
One  assumes  that  the  wind  speed  is  zero  and  that  the  relative 
humidity  is  100%  at  the  surface.  The  air-surface  differences  in 
these  bulk  parameters  are  used  to  calculate  surface  layer  fluxes 
and  stability.  The  measurements  must  be  made  in  the  "surface 
layer",  which  is  a  layer  approximately  30m  thick,  where  the  heat 
flux  and  friction  velocity  can  be  assumed  constant. 

The  Monin-Obukhov  length  is  calculated  from 

L  *  (T/kg) (U*3/q) ,  (1) 

where  T  is  the  absolute  temperature,  k  is  Von  Kantian* s  constant 
(0.35),  g  the  acceleration  due  to  gravity,  and  U*  the  friction 
velocity  (scaling  velocity).  The  temperature  flux,  Q,  is  related 
to  the  virtual  heat  flux,  H,  by 

H  *  pCpQ,  (2) 

where  p  is  the  air  density  and  Cp  is  the  specific  heat  at 
constant  pressure.  Note  that  for  upward  heat  flux  (unstable  air) 
L  is  negative  and  is  positive  for  downward  heat  flux  (stable  air) 
We  will  work  with  l/L  in  what  follows,  which  has  the  value  zero 
for  neutral  conditions  where  the  heat  flux  is  zero.  The 
temperature  flux  is  directly  related  to  the  virtual  potential 
temperature  scaling  parameter,  9y*,  by 

Q  *  ^  2 j 

■  u*(e*  +  6.1  x  10  "4  Tq* ) , 
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where  q*  is  the  water  vapor  mixing  ratio  scaling  parameter 
(gm/kg).  The  potential  temperature  is  9  *  T  +  0.G098Z.  We  see 
that  L  can  be  calculated  from  the  surface  layer  fluxes,  or  the 
scaling  parameters,  since  they  are  directly  related. 

Briefly,  the  method  we  use  is  to  determine  the  scaling 
parameters  directly  from  the  measured  mean  values  using  stability 
corrected  drag  coefficients.  Note  that  this  leads  immediately  to 
an  iterative  calculation  since  the  stability  must  be  known  to 
determine  the  proper  drag  coefficient,  stability  (L)  is  determined 
from  the  scaling  parameters,  and  the  scaling  parameters  depend  on 
the  drag  coefficient.  The  procedure  is  fairly  simple  and  is 
described  in  detail  in  Appendix  A. 

B  .  Rich ar d-aon  '  s  N umber 

We  do  not  use  the  Richardson's  number  in  what  follows,  but  since 
it  is  often  used,  we  present  it  here  for  the  sake  of  completeness. 
Like  the  Monin-Obukhov  length  it  can  be  determined  from  either 
bulk  or  gradient  measurements.  The  gradient  Richardson's  number, 
determined  from  a  two  level  measurement  is 

Ri  *  Z(g/T)C(92  -  91)2/(U2  *  Ul)]ln(Z2/Z1),  (4) 

where  the  subscripts  1  and  2  refer  to  specific  heights  and  Z 
is  the  geometric  mean  height.  The  bulk  Richardson’s  number  is 
determined  from 

B  *  Z ( g/T) (9  -  0O)/U2,  (5) 
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where  0O  is  the  surface  value.  Ri  and  L  can  be  approximately 
related  to  each  other  by 

l/L  =  Ri/Z  (unstable)  (6) 

and 

l/L  *  Ri / Z ( 1  -  SRi)  (stable)  (7) 

The  value  of  8  is  taken  to  be  7,  although  this  parameter  is  not  well 
known. 


C.  Pasquill  Class 

The  Pasquill  classification  is  a  simplified  scheme  which 
requires  a  determination  of  wind  speed  but  not  temperature. 

A  heat  flux  category  is  determined  from  an  estimation  of  the  solar 
insolation  at  the  surface.  The  scheme  has  been  developed  for 
overland  use  and  is  not  appropriate  for  the  overwater  regime 
because  of  the  vastly  different  thermal  properties  of  water. 
Stability  is  divided  into  five  classes,  A-F,  determined  as 
follows: 


U (m/sec) 
<2 
2 
4 
6 

>6 


Daytime  Insolation 


Strong 

A 

A-B 

B 

C 

C 


Moderate 

A-B 

B 

B-C 

C-D 

D 


Small 

B 

C 

C 

D 

D 


Nightime  Cloud  Cover 
>4/8  <3/8 


E 

D 

D 

D 


F 

E 

D 

D 
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In  this  scheme  A-C  are  unstable,  D  neutral,  and  E-F  stable 


conditions.  This  classification  scheme  is  imprecise  but,  due  to 
its  simplicity,  is  in  widespread  use.  One  of  our  purposes  here  is 
to  deduce  a  similar  scheme  for  use  in  the  overwater  regime,  and 
which  will  be  more  precise  than  the  Pasquill  categories. 

D.  EPA  Temperature  Gradient  Classification 

The  Environmental  Protection  Agency  has  guidelines^)  for 
relating  the  Pasquill  classifications  to  the  wind  direction 
standard  deviation,  a @,  and  also  to  the  temperature  gradient 
in  the  surface  layer.  This  scheme  considers  only  the  temperature 
gradient;  no  determination  of  the  wind  is  used.  Categorization  is 
as  follows: 


AT/AZ( °C/100m) 

Pasquill  Class 

Ja(deq) 

<-1.9 

A 

>2  2.5 

-1.9  -v  -1. 7 

B 

22.5  17.5 

-1.7  -  -1.5 

C 

17.5  -  12.5 

-1.5  -  -0,5 

D 

12.5  *  7.5 

-0.5  -  1.5 

E 

7.5  *  3.8 

1.5  -  4.0 

F 

<3.8 

>4.0 

G 

It  is  important  to  realize  that  this  scheme  was  also  developed  by 
examination  of  overland  data;  it  is  not  expected  to  be  valid  over 


water . 


V.  Modification  of  Pasquill  Classification  for  Overwater  Use 

As  has  been  stated  above,  the  Pasquill  classification  scheme 
was  derived  for  overland  use.  The  overwater  regime  is  different 
because  of  the  vastly  different  response  of  water  to  incoming 
solar  radiation  and  the  variation  of  roughness  length  with  wind 
speed.  Because  of  the  widespread  use  of  this  scheme  we  here 
modify  it  for  overwater  use;  this  should  not  be  interpreted  as  a 
recommendation  that  this  scheme  remain  in  common  use . 

Golder^4)  has  related  the  Pasquill  class  boundaries  to 
the  roughness  length  and  the  Monin-Obukhov  length.  The  results 
are  shown  in  Figure  3.  The  dependence  of  the  roughness  length  on 
wind  speed  has  received  much  attention.  We  have  chosen  the 
parameterization  of  Kondo^),  which  is  detailed  in  Appendix  A. 
Figure  3  has  both  wind  speed  and  roughness  length  scales;  the  wind 
speed  scale  was  derived  from  this  parameterization. 

Wind  speed,  temperature,  and  relative  humidity  (or  dew  point) 
are  the  meteorological  parameters  which  jure  commonly  measured,  not 
1/L.  Of  course,  1/L  can  be  determined  from  the  mean  parameters  by 
the  method  described  in  Appendix  A.  We  have  calculated  1/L  as  a 
function  of  wind  speed,  air-sea  temperature  difference  (TA-TS), 
and  relative  humidities  of  50%,  80%,  and  95%.  The  results  are 
shown  in  Figures  4.  Also  plotted  on  these  figures  are  the 
boundaries  of  the  Pasquill  classes,  as  determined  from  Figure  3. 
Note  that  classes  A  and  F  do  not  appear  on  Figures  4  because  these 
classes  are  outside  the  range  of  parameters  that  can  be  expected 
to  occur  over  water  (particularly,  very  large  air-sea  temperature 
difference) . 

Finally,  using  Figures  4,  we  plot  the  Pasquill  class 
boundaries  as  functions  of  wind  speed  and  TA-Tg  in  Figure  5. 


18 


There  are  three  pLots,  one  each  for  50%,  80%,  and  95%  relative 
humidity.  The  humidity  is  not  a  major  factor;  it  is  most 
important  near  neutral  (small  TA-Tg)  where  high  humidity 
shifts  the  curves  toward  greater  instability. 

Note  that  all  of  the  above  calculations  have  been  made  for  a 
height  of  20m,  for  the  measured  mean  parameters. 

Because  of  the  wide  range  of  parameters  applicable  for  each 
class,  particularly  D,  we  have  arbitrarily  divided  each  class  into 
three  subclasses,  each  subclass  representing  one  third  of  the 
original  class.  For  example,  D+,  D,  and  D-  where  the  +/-  refer  to 
the  stable/unstable  side  of  the  original  class. 
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u  (m/sec) 


(oas/ui)  n 


u  (m/sec) 


/sec)  u  I  m/sec) 


The  ultimate  test  of  any  classification  scheme  is  how  well  it 
can  predict  the  behavior  of  a  plume.  Such  a  comparison  for  the 
data  set  gathered  in  this  experiment  is  being  made  by  AI.  An 
equivalent  test  can  be  made  by  comparing  to  the  variability  of  the 
wind  direction,  which,  of  course,  is  the  major  cause  of  cross  wind 
spreading  of  a  plume.  This  is  the  approach  taken  here.  Short 
term  (high  frequency)  fluctuations  will  cause  spreading  of  a  plume 
about  its  center  line  while  long  term  variations  will  cause  plume 
meander.  One-hour  average  measurements  of  the  plume  concentration 
profile  and  one-hour  averages  of  the  wind  direction  standard 
deviation  will  both  include  short  and  long  term  variability.  In 
this  report,  we  will  make  only  comparisons  for  long  term  averages 
since  the  purpose  here  is  to  parameterize  models  that  compute 
one-hour  average  impact.  The  data  set  used  (shown  in  section  III) 
includes  periods  when  tracer  gas  was  not  released.  Measured 
values  of  Oq  at  a  single  point  can  be  used  to  calculate  the 
plume  profile  only  if  horizontal  homogeneity  exists. 

We  have  compared  the  measured  values  of  3  0  with 

A.  air-sea  temperature  difference  and  temperature  gradient 

B.  10/L 

C.  modified  Pasquill  class 

The  air-sea  temperature  difference  and  temperature  gradient  are 
examined  because  of  the  EPA  quidelines.  We  do  not  expect  them  to 
be  of  much  use  because  stability  cannot  be  determined  from 
temperature  information  alone. 


0.4  -0.2 


0  0.2  0.4 

T20-T7  PC) 


A.  Air-Sea  Temperature  Difference  and  Temperature  Gradient 


On  the  ship  the  air  temperature  at  7m  and  20m  above  sea 
level,  and  the  sea  surface  temperature,  Tg,  were  measured 
continuously  and  averaged  every  half  hour.  A  description  of  the 
measurements  and  equipment  used  is  given  in  Appendix  B.  For 
air-sea  temperature  difference  we  used  (T2Q-T5)  and  for 
temperature  gradient  we  used  (T'20_T7)  directly,  not  dividing 
by  the  height  difference.  versus  (T2o~T7^  and 

(t2o-tS>  are  shown  in  Figures  6  and  7,  respectively. 

There  is  little  correlation  between  aQ  and  these 
parameters  so,  as  expected,  it  would  be  very  difficult  to 
determine  a@  from  temperature  measurements  alone.  For  the 
sake  of  completeness  we  did  a  linear  regression  for  both  sets  of 
data  for  negative  temperature  gradient.  There  is  not  enough  data 
for  negative  temperature  gradient.  There  is  not  enough  data  to 
warrant  the  calculation  for  positive  gradient.  The  correlation 
coefficients  of  the  linear  regressions  were  only  0.14  for  the 
gradient  and  0.18  for  the  air-sea  difference.  It  is  interesting 
to  note  that  (T20“Ts)  values  range  from  -2  to  +2  C*  whereas 
the  range  of  (^20~^7^  is  a  factor  of  ten  smaller.  This 
emphasizes  the  difficulty  in  measuring  temperature  differences  in 
the  air  over  the  ocean.  The  differences  are  normally  too  small  to 
yield  reliable  results.  It  is  much  better  to  use  the  air-sea 
temperature  difference  when  stability  is  being  determined. 

B.  10/L 

We  have  compared  hourly  averages  of  the  one  minute,  ten 
minute  and  one  hour  averages  of  aQ  with  hourly  averages  of 
10/L.  The  results  are  shown  in  Figure  8. 


For  unstabLe  conditions  somewhat  better  results  were  obtained  when 


Jq  versus  ln(-10/L)  was  plotted  and  this  is  shown  in  Figure 
9.  For  both  figures  the  solid  Lines  are  linear  regression  fits  to 
the  data.  The  open  circles  in  Figure  9  are  the  averages  for  all 
data  in  the  bins  -10/L  <0.1,  0.1  _<-10/L  <  0.2,  0.2  <_  -  10/'L  <  0.5, 
0.5  <_  -10/L  <  1.0,  1.0  -10/L.  The  linear  regression  results 

are 


Plot 


intercept 


aQ  vs  — 1 0/L ,  unstable  5.4 
Og  vs  10/L,  stable  9.7 
CJg  vs  ln(-10/L),  unstable  9.2 


s lope  correlation 

2.7  0.45 

26  0.47 

1.9  0.57 


The  correlations  are  not  good,  but  they  are  not  surprising  when 
one  considers  the  wide  range  of  conditions  encountered  and  the 
region  in  which  the  measurements  were  made.  We  were  not  dealing 
with  a  spatially  homgeneous  air  mass  which  had  a  long  fetch  over 
the  open  ocean.  As  described  above,  the  area  has  a  complex 
pattern  of  eddies  and  channeling  of  the  wind.  With  the  limited 
data  set  we  have  available  it  is  not  possible  to  bin  the  data  into 
categories  which  delineate  the  local  flow  conditions,  which  might 
reduce  some  of  the  observed  scatter. 

The  fact  that  ag  shows  a  better  correlation  with 
ln(-10/L)  than  with  -10/L  is  rather  peculiar,  and,  in  fact,  we 
include  Figure  9  only  because  it  better  shows  that  og  is 
correlated  with  the  Monin-Obukhov  length.  It  would  be 
unreasonable  to  use  a  logarithmic  dependence  on  stability,  unless 


-10/L  (m-1) 


the  functions  were  truncated,  since  jg  would  become  negative 
as  neutrality  is  approached. 

We  have  used  the  above  results  to  determine  average  relations 
for  0q  as  a  function  of  10/L.  They  are  ( in  degrees, 

L  in  meter  s ) : 
stable  conditions 

Jq  =  12  (10/L  >  0) 

unstable  conditions 

0Q  =  5.4  +  2. 7( -10/L)  (10/L  <  0) 

or 

JQ  =  10.2  +  2.4  ln( -10/L)  (10/L  <  -0.05) 

=  3.0  (-0.05  <_  10/L  <  0) 

This  combination  of  relations  is  rather  peculiar  since  they  do  not 
agree  at  10/L  *  0.  The  linear  regression  fit  to  10/L  does  not 
appear  to  do  a  good  job  near  neutrality  for  unstable  conditions 
and  is  probably  meaningless  for  the  small  data  set  for  stable 
conditions  data  set.  The  truncated  logarithmic  fit,  and  the 
constant  value  for  stable  conditions,  best  approximate  the  data  we 
have  obtained.  Note  in  Figure  9  that  the  open  circles,  which  cure 
the  average  values  for  the  ln(-10/L)  bins  listed  above,  are  fairly 

4 

well  approximated  by  the  linear  regression  fit  (solid  line).  One 
important  note  must  be  made:  when  using  Equations  8  and  10  the 
calculated  values  of  can  be  in  error  by  as  much  as  a 
factor  of  ten  and  errors  of  a  factor  of  two  will  be  common. 

We  have  attempted  to  reduce  the  scatter  in  the  results  by 
plotting  Oq  versus  wind  direction  and  10/L  to  see  if  the 
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Local  air  mass  trajectory  had  an  effect.  The  results  are  shown  in 
Figure  10,  where  different  symbols  are  used  for  four  ~q 
bins.  The  results  are  not  helpful;  wind  direction  does  not  appear 
to  play  an  important  role. 

The  most  striking  aspect  of  these  results  is  the  large 
for  slightly  stable  conditions.  Previous  experiments, 
over  both  land  and  water,  have  shown  a  continuous  increase  of 
3q  (or  Oy)  as  atmospheric  stability  goes  from  stable  to 
unstable  (6,7,8) ^  we  believe  the  results  of  this  study  are 
location  specific.  Stable  conditions  over  California  coastal 
waters  are  unusual,  and  normally  associated  with  offshore  flow  or 
where  the  air  mass  has  not  equilibrated  after  being  transported 
across  a  change  in  sea-surface  temperature.  In  the  Ventura  area, 
local  eddies  or  flow  over  the  channel  islands  could  produce  a 
stable  surface  layer.  In  both  cases,  variations  in  the  wind 
direction  could  be  large  due  to  orographic  effects.  Such 
variations  would  presist  over  long  distances  due  to  the  stability 
of  the  atmosphere.  Thus,  the  stable  results  reported  here  for 
3q  should  not  be  thought  to  apply  to  stable  onshore  flow  in 
open  coastal  areas. 

C .  Overwater  Modified  Pasguill 

We  have  determined  one-hour  averaged  values  of  the  modified 
classes  from  the  measurements  made  on  the  RV/Acania,  using  the 
method  described  in  Section  V.  The  measured  values  of  3q 
and  the  modified  classes  for  these  time  periods  are  shown  in 
Figure  11.  Even  though  there  are  48  hours  of  data,  there  is 
insufficient  data  to  yield  a  good  versus  class 


relationship  because  so  much  of  the  data  are  near  neutral.  We 
cannot  verify  that  the  classes  should  be  divided  into  subclasses, 
though  it  does  appear  to  be  valid  for  class  D.  The  results 
suggest  the  following  relation  to  obtain  the  classes  from 
a@,  or  vice-versa: 


Pasquill  Class 


10  -  14 


6-10 


D  (unstable) 


2-6 


D  (stable) 


7-14 


•  Unstable 

*  Stable 


APPENDIX  A 


Determination  of  Stability  and  Scaling  Parameters 


We  develop  here  the  method  used  to  determine  stability  from 
measured  mean  meteorological  parameters.  The  method  also  yields 
the  surface  layer  scaling  parameters,  which  are  directly  related 
to  the  fluxes.  We  also  show  how  the  convective  mixing  velocity 
can  be  determined.  This  velocity  can  be  used  to  determine  the 
rate  at  which  pollutants  are  mixed  through  the  full  boundary 
layer . 

In  what  follows  we  use  the  notation: 

T  temperature  (K) 

Z  height  (m) 

U  mean  wind  speed  (m/sec) 

q  water  vapor  mixing  ratio  (g/kg) 

9  =  T  +  . 0098Z  potential  temperature 

9y  —  9  +  .61Tq  virtual  potential  temperature 

g  acceleration  due  to  gravity 

k  von  Karman's  constant  (0.35) 


U* 

q* 

T* 

w* 

zo 

L 


friction  velocity  (scaling  velocity) 
mixing  ratio  scaling  parameter 
scaling  temperature 
convective  mixing  velocity 
roughness  length 
Monin-Obukhov  length 


X 

ax 

X 

C 

«5) 

♦(5) 


inversion  height 

atmospheric  parameter  (T,  q  or  U ) 
turbulent  diffusivity  ratio  ( ctjj  =  1) 
drag  coefficient 
Z/L 

gradient  stability  function 
profile  stability  function 


We  describe  the  stability  with  the  Monin-Obukhov  length 


L 


T  U*2 
_ _ • 

©v* 


(A-l) 


The  scaling  parameters  are  related  to  the  gradients  by 


dX  _ 
dZ 


a  kZ 
X 


*X  (5)- 


(A-2) 


Current  best  evidence  shows  that  transport  of  the  scalars,  heat 
and  water  vapor,  obey  the  same  relationship  to  their  gradients. 
Thus,  T  *  4>q  and  aT  *  ciq  *  1.35.  The  stability 
correction  functions  are: 


♦  t(5)  -  (1  -  9C)'1/2  S<  0  u<  C)  =  (1  -  15C)"1/4 


(A- 3) 


-  (1  +  6.4  )  5  >  0  =  (1  +  4.70 
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BuLk  Method 


To  obtain  the  stability,  the  scaling  wind  speed  and 
temperature  are  needed  for  use  in  Equation  Al.  They  are  obtained 
from  bulk  measurements  using  an  integrated  form  of  Equation  A2. 

We  integrate  Equation  A2  from  the  surface  to  a  reference  height  Z, 
usually  10m- 


xz  -  xs 


- -/ 

ak  J 


Z  4>x (? ) 


dZ 


(A-4) 


X*  Z 

=  —  [in  - 

ak  Z_ 


+'x(5  5  1  / 


where  Xs  is  the  surface  value.  We  have  dropped  the  subscript  on 
a  for  convenience. 

For  analysis  of  data,  it  is  most  convenient  to  solve  Equation 
A4  for  X*  and  rewrite  in  terms  of  the  drag  coefficient.  The 
neutral  stability  drag  coefficient  is  given  by 


*4  =  3 

XN  110  z/zoX 

and  corrected  for  stability  by 

X"  NX^1  ■*  Hi*  )CNx1/2/ak]-1 
Thus,  the  scaling  parameter  is  given  by 


(A- 5) 


(A-6) 
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(A-7) 


X*  =  C*X(X2  -  Xg) 


As  was  stated  above,  the  key  to  the  analysis  is  obtaining  the 
stability.  From  Equation  1  we  write  £  as 


5  =  k  g  Z  8  +  0.61  x  10~4  Tq* 
T  U*2 


(A-8) 


For  the  humidity  correction  term,  we  use  T  =*  15*C  and  approximate 
0.61  T  3  0.18  (q  in  gm/kg).  Rewriting  in  drag  coefficient  form 
gives,  for  neutral  stability, 


r 

"o 


kg  Z  TtN  (T  -  T«)  t  0.18  (q  -  qa ) . 
T  ^JN  U2 


(A-9) 


We  assume  zero  wind  speed  at  the  surface.  Using  Equations  A6,  A8, 
and  A9  gives 

5  =  5q  t  1  ^u(C)Cujj/k]  (A-10) 

1  -  V  -p  (  5 ) Ctn/  °k 

Obviously,  in  order  to  determine  £  the  stability  functions  must 
be  known  and  they  are  determined  from  £ .  Thus,  an  iterative 
procedure  must  be  used.  The  steps  are  as  follows: 

1.  Calculate  q  and  qg  from  the  measured  Tg,  T,  and  dew 

point  temperature  for  Z  «  10  m  and  assuming  100%  humidity 
at  the  surface. 
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r 


2.  Use  =1.3  x  10~3  an<j[  from  the  table 

below  and  the  measurements  to  obtain  £Q. 

3.  Calculate  'i  ^  and  by 

4.  Calculate  £  from  Equation  A10. 

5.  Iterate  steps  3  and  4  until  the  desired  accuracy  is 
obtained,  giving  £ . 

6.  L  has  been  determined  and  U*,  T*,  and  q*  are 
obtained  directly  from  Equation  A7. 

For  Z  =  10  m  the  wind  drag  coefficient  is  found  from  5 


U  (m/sec)  ^iw  x  103 


0.3 

-  2.2 

1.08 

c 

1 

h-* 

cn 

2.2 

-  5.0 

0.77 

+  0.086U 

5  - 

8 

0.87 

+  0.067U 

8  - 

25 

1.2 

+  0.025U 

The  profile  stability  functions  are 
■pT(  £<  0)  =  2in  (i~)  for  x  *  (1  -  95)1/2 

(  5  >  0)  =  -6.55 

2 

P  U  (  5  <  0)  *  2  £  n  ( -~2~)  +  £  n(  — ^  - — )  —  2  tan  ^x  +  j- 
for  x  «  (1  -  15  5)1/4 


(A-ll) 


(A-12) 


*D<  5  >  0)  *  -4.75 
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Surface  Layer  Fluxes 


The  scaling  parameters  described  above  apply  to,  and  specify 
the  state  of  the  surface  layer.  The  surface  layer  momentum,  heat, 
and  water  vapor  fluxes  can  be  determined  from  scaling  parameters 
by 

Fm  «  PU*2 

Fh  *  pCpU*  9*  (A- 13) 

Fq  *  pU*  Q* 

The  F-fj  in  Equation  A13  is  the  sensible  heat  flux.  The  depth 
over  which  the  surface  scaling  applies  is  approximately  |  L  |. 

Boundary  Layer  Scaling 

Above  the  surface  layer  (ignoring  the  transition  region  to 
the  well  mixed  layer) ,  surface  scaling  no  longer  applies  and  new 
scaling  length,  velocity,  and  temperature  are  needed.  These 
are:  ^ 

Zi 

W*  =■  C  ( g/T )  Qq  Zili/3  (A-14) 

H*  »  Q0/w* 

We  assume  that  the  boundary  layer  depth  is  defined  by  the 
inversion  height,  Z^.  Q0  is  the  surface  virtual  heat  flux. 

w*  is  the  scaling  velocity  in  the  well  mixed  layer,  and  we 
assume  that  this  velocity  is  directly  related  to  the  boundary 
layer  convective  mixing  velocity.  Previous  SFg  tracer 
experiments  performed  by  California  Institute  of  Technology  in 
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cooperation  with  NPS  have  shown  that  w*  computed  by  the  above 
method  closely  predicts  the  tracer  experiment  results.  Laboratory 
experiments  (10)  have  also  shown  that  the  vertical  transport 
of  particulate  matter  is  given  by  w* . 

Equation  A-14  for  w*  is  not  correct  for  the  cloud  topped 
boundary  layer.  In  that  case  the  presence  of  the  cloud  modifies 
the  heat  flux  due  to  release  of  latent  heat  upon  condensation  and 
absorption  and  emission  of  radiation.  The  term  Q0  in  Equation 
A-14  must  be  replaced  with  2.5  <Q>,  the  average  heat  flux  in  the 
boundary  layer.  Use  of  the  equation  as  it  stands  will  typically 
lead  to  underestimates  of  about  30%  {it  is  this  small  due  to  the 
one- third  power  dependence  of  w*  on  Q) . 
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APPENDIX  B 


Shipboard  Meteorological  Equipment  and  Calibration  Procedures 

The  RV/Acania  has  two  meteorological  stations  at  heights  of 
7.0  and  20.5  meters  above  mean  sea  level  and  a  wide  range  of 
auxiliary  instrumentation.  Sensors  are  mounted  on  masts  in 
foreward  locations  so  that  ship  disturbance  of  the  airflow  has  as 
little  influence  as  possible.  We  will  not  describe  all  of  the 
shipboard  equipment,  only  that  which  was  used  to  gather  the  data 
for  this  study.  A  complete  description  appears  elsewhere ( 11 ) . 

The  meteorological  stations  have  equipment  to  measure  the 
following  parameters: 

relative  wind  speed 
relative  wind  direction 
air  temperature 
dew  point 

Wind  direction  was  measured  only  at  the  top  level.  The  true  wind 
speed  and  direction  were  calculated  from  the  relative  wind  speed 
and  direction  and  the  ship's  speed  and  heading.  Auxiliary 
measurements  which  were  made  included: 

sea  surface  temperature 
inversion  height 
ship  roll 
ship  heading 
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Individual  Instrumentation 


Wind  speed  and  direction;  Meteorology  Research,  Inc.,  1022 
sensors  with  1001  transmuter.  The  wind  direction  sensor  has  dual 
potentiometers?  the  transmuter  has  automatic  switching,  giving  a 
0-540°  continuous  readout.  The  wind  speed  sensors  are  cups  with 
a  0.25  m/sec  threshold,  60  m/sec  maximum  speed,  and  1%  accuracy. 
Dew  Point;  General  Eastern,  1200AP,  with  1211P  sensor.  The 
system  measures  the  dew  point  by  the  cooled  mirror  technique.  The 
change  in  the  reflectivity  of  a  cooled  mirror  upon  condensation  is 
detected  optically.  The  optical  detector  is  in  a  servo  circuit 
which  controls  the  mirror  temperature,  holding  it  at  the  dew 
point.  A  platinum  resistor,  bonded  to  the  mirror,  is  the 
temperature  detector.  The  sensor  is  mounted  in  an  aspirated 
radiation  shield  suppled  by  the  manufacturer.  The  system  has  a 
range  of  approximately  70 °C  below  ambient  temperature,  which  is 
governed  by  the  capacity  of  the  thermoelectric  cooler.  Operating 
temperature  limits  are  100°C  to  -59°C.  Manufacturer  stated 
accuracy  is  +  0.22°C. 

The  platimum  thermometer  comes  equipped  with  a  three  wire 
hook  up  for  the  readout  that  is  supplied  with  the  instrument.  The 
readout  converts  sensor  resistance  (three-wire)  to  a  voltage 
proportional  to  temperature.  Because  of  drifts  in  the  voltage 
conversion  circuitry,  we  read  the  sensor  resistance  with  a 
computer  controlled  4-wire  Ohmeter  and  convert  to  temperature  in 
computer  software.  Converting  the  sensor  to  a  4-wire  system  has 
eliminated  lead  resistance  influence. 


Air  temperature:  Rosemount,  100  Ohm  platinum  resistors,  4-wire 
configuration.  The  resistance  is  measured  with  the  same  computer 
controlled  meter  used  for  the  dewpoint  measurements.  The  resistor 
is  mounted  in  a  1/4  in  stainless  steel  tube.  This  tube  fits 
loosely  in  a  3/8  inch  tube  topped  by  a  small  metal  enclosure.  The 
enclosure  has  an  environmentally  tight  screw  cap  so  it  may  be 
opened  for  access  to  the  resistor  leads.  The  outer  tube  and 
housing  are  enclosed  in  an  R.M.  Young,  Gill  aspirator  for 
radiation  shielding.  The  aspirator  consists  of  a  silvered  glass 
dewar  which  is  isolated  from  the  heat  of  the  aspiration  motor  by  a 
long  tube.  To  improve  the  performance  of  the  aspirator  we  have 
mounted  a  silvered  mylar  radiation  shield  at  the  inlet.  This  has 
been  found  to  be  especially  necessary  on  shipboard  where  some 
sensors  are  over  the  water  and  some  over  the  warm  deck  of  the 
ship. 

In  order  to  eliminate  contact  e.m.f.  problems  the  4  wires 
leads  are  continuous  all  the  way  from  the  resistor  to  the  Ohmeter . 
Solder  connections  are  made  in  the  aspirator  housing.  All  solder 
points  are  at  the  same  temperature  so  that  thermal  e.m.f.  is  not  a 
factor.  The  accuracy  of  the  system  is  +0.2*C,  the  uncertainty 
being  due  to  imperfect  aspiration.  In  order  to  achieve  this 
accuracy  with  the  100  Ohm  platinum  resistor  the  Ohmeter  must  have 
5  significant  figure  accuracy. 

Sea-Surface  Temperature  :  The  sensor  is  a  100  Ohm  platinum 
resistor  identical  to  that  used  for  air  temperature.  The  sensor 
is  enclosed  in  a  long  streamlined  float  that  can  be  towed  at 
normal  ship  speed  or  will  remain  near  the  surface  when  the  ship 
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is  stopped.  The  float  is  a  12  ft  long,  3/4  inch  internal  diameter, 
tygon  tubing.  One  end  of  the  tubing  is  sealed  with  a  brass  plug 
and  the  sensor  is  mounted  with  good  thermal  contact  in  this  plug. 
The  plug  extends  1  inch  beyond  the  end  of  the  tubing  so  that  it 
makes  good  thermal  contact  with  the  water.  The  other  end,  through 
which  the  lead  passes,  is  sealed  by  silicone  cement.  The  sensor 
is  suspended  from  the  side  of  the  ship  on  a  guyed  pole  so  that  it 
is  outside  the  ship's  influence.  The  sensor  samples  approximately 
the  top  12  inches  of  the  water  due  to  its  bobbing  with  ship  motion 
Accuracy  is  +  0.1*C. 

Inversion  Height:  Aerovironment  300,  Acoustic  Radar  with  100  Watt 
driver.  The  device  has  a  range  of  100  to  1000  m,  but  it  is 
difficult  to  discern  an  inversion  above  700  m.  Shipboard  noise 
limits  the  useful  range  to  about  500  m  when  the  ship  is  underway. 
The  enclosure  used  for  the  antenna  is  a  pentagon  with  sides  4  ft 
wide  and  8  ft  high.  For  shipboard,  the  base  of  the  enclosure  must 
be  soundproofed.  The  recently  available  toothed  enclosure 
extensions,  Thenadners,  also  reduce  the  noise  considerably.  The 
enclosure  is  normally  located  near  the  stern  of  the  ship  to  reduce 
wind  noise,  and  in  a  location  where  stack  noise  is  minimal. 

Ship  Heading;  On  the  RV/Acania  a  Librascope  7070  analog  to 
digital  converter  has  been  connected  to  the  gyroscope  used  for 
navigation.  The  converter  is  a  mechanical  device  that  senses  the 
readout  angular  orientation.  A  buffer  interfaces  the  output  of 
the  converter  to  a  16  bit  interface  to  the  computer. 
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Ship  Roll  and  Roll  Rate:  A  mechanical  pendulum  geared  to  a  poten¬ 
tiometer  placed  on,  and  pivoting  about,  the  roll  axis  of  the  ship, 
provides  an  analog  signal  proportional  to  the  ship's  roll  angle. 

A  simple  RC  network  is  used  to  differentiate  this  signal  to  yield 
an  output  related  to  roll  rate.  Errors  introduced  by  vertical  and 
lateral  accelerations  of  the  ship  are  minimized  by  placing  the 
pendulum  on  the  roll  axis  and  midway  along  the  length  of  the 
ship. 

Data  Acquisition  Equipment 

Three  basic  methods  of  data  acquisition  and  shipboard  analysis 
are  used:  1)  Computer  controlled  acquisition  of  all  signals,  2) 
Analog  tape  recording  and  strip  charts,  and  3)  Spectral  analysis. 
Spectral  Analysis;  Nicolet  Scientific  Corp.,  440B  Spectrum  analy¬ 
zer.  We  have  found  that  spectral  analysis  is  a  necessity  on 
shipboard,  particularly  to  monitor  fluctuation  signals.  A  ship  is 
an  electrically  noisy  environment,  and  the  noise  is  not  constant. 
Spectral  analysis  can  help  solve  noise  problems.  Frequent 
monitoring  of  signals  to  verify  the  quality  of  the  data  is 
advisable. 

Analog  recording:  Honeywell  5600E  tape  recorder.  This  unit  has 
14  analog  channels  and  automatic  speed  jontrol.  We  have  built  a 
selector  so  that  all  14  input  or  output  channels  can  be  monitored 
on  a  oscilloscope  and  by  spectral  analysis.  The  recorder  is  a 
back  up  device,  all  analog  signals  being  recorded  in  case  there  is 
a  failure  in  any  other  portion  of  the  data  acquisition  system. 


48 


Strip  chart  recordings  are  made  of  the  relative  wind  direc¬ 
tion,  and  any  three  other  parameters,  as  desired.  The  strip  chart 
allows  post-experiment  analysis  of  the  data  to  be  done  to  pick  out 
times  when  external  influence  (such  as  RFI)  could  invalidate  the 
data.  The  strip  charts  allow  these  events  to  be  detected 
(recognizing  that  qualified  personnel  may  not  be  available  for 
continuous  24-hour  monitoring) . 

Computer  controlled  acquisition;  Hewlett  Packard  (HP)  9825  or 
9835  desktop  computer  and  HP  3052A  data  acquisition  system.  The 
data  acquisition  system  includes  a  59309A  digital  clock,  3455A 
digital  voltmeter,  and  3495A  scanner.  Data  output  is  on  the 
computer  internal  tape  and  a  9871A  or  a  2631  printer. 

The  scanner  has  40  channels  of  low  thermal  emf  relays  with  a 
switching  time  of  less  than  10  msec  and  a  contact  differential  emf 
of  less  than  2  yVolt.  The  relays  are  three  terminal, 
accommodating  a  hi,  lo,  and  ground  or  guard.  This  is  desirable 
since  the  guards  for  the  various  signal  channels  can  be  kept 
independent,  eliminating  ground  loops. 

The  voltmeter  has  six  significant  figures  resolution,  which  is 
needed  for  the  platinum  resistance  thermometer  measurements. 
Accuracy  with  the  high  resolution  off  is  0.001%. 

The  interface  between  the  computer,  data  acquisition  system, 
and  printer  is  the  HPIB  (IEE-488-1975 ) .  In  addition,  the  computer 
is  interfaced  with  the  buffer  for  the  LORAN-C  navigation  and  the 
ship  heading  indicator  by  a  16  bit  interface.  The  computer 
alternately  operates  the  scanner,  interrogates  the  voltmeter,  and 


49 


interrogates  the  buffer.  Ail  scanner  channels  are  activated  and 
the  data  acquired  between  individual  readings  of  the  heading. 

Normal  shipboard  procedure  is  to  produce  half-hour  averages  of 
all  signals.  The  exact  number  of  samples  of  each  channel  depends 
on  the  program  being  used,  e.g.  some  programs  do  many  more 
calculations  and  intermediate  averages.  The  number  of  samples 
ranges  from  800  to  1500.  The  computation  plus  print  time  uses  one 
to  two  minutes  of  each  half-hour  period.  During  the  half-hour 
cycle  all  signals  are  acquired,  averaged,  computations  made,  and 
all  averaged  voltages  and  computed  parameters  are  stored  on  the 
tape  and  printed  out.  Included  in  the  storage  and  printout  are 
all  system  calibration  parameters  so  that  all  calculations  can  be 
redone  later,  if  needed. 

Calibration  Procedures 

The  results  of  the  laboratory  calibrations  of  meteorological 
sensors  for  use  in  field  measurements  can  be  misleading.  For 
example,  we  can  calibrate  the  platinum  resistance  temperature 
sensors  to  +  0.001°C,  or  can  calibrate  them  with  the  leads  used 
in  the  field  to  +  0.01°C.  This  calibration  is  repeatable  and 
represents  the  basic  capabilities  of  the  calibration  facility. 

When  these  sensors  are  used  in  the  field,  with  the  required 
aspiration,  recognizing  that  at  times  the  aspirator  will  be  shaded 
and  at  times  in  direct  sunlight,  measurement  accuracy  of  +  0.1 “C 
is  doing  well.  Measuring  wind  speed  is  another  problem  because  of 
ship  influence.  No  matter  how  accurately  the  cups  are  calibrated 
using  them  on,  or  near,  the  RV/Acania  reduces  the  accuracy  to 
+0.5  m/sec.  The  results  on  larger  ships  are  worse.  The 
attainable  accuracies  in  the  field  are  sufficient  for  this  work. 
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especially  when  one  considers  the  use  of  the  results. 
Meteorological  parameters  measured  at  a  point  are  being  used  to 
predict  plume  behavior  over  a  spatially  large  region,  near  the 
coast,  where  spatial  inhomogeneities  are  great.  Great  accuracy  is 
not  required  but,  of  course,  the  sensors  should  be  moderately  well 
calibrated  so  that  systematic  errors  are  not  introduced.  As  will 
be  seen  in  the  following  description  of  calibration  procedures,  we 
are  well  within  that  goal. 

Air,  Sea  and  Dewpoint  Temperatures:  All  temperature  sensors  are 
calibrated  together  in  an  insulated  chamber.  The  air  and  sea 
temperature  probes  and  a  HP2801A  quartz  thermometer  standard  are 
mounted  in  an  aluminum  cylinder  which  effectively  eliminates  the 
temperature  differences  between  these  sensors.  The  differences 
are  less  than  0.01 °C.  The  dewpoint  thermometers  cannot  be  mounted 
in  the  same  block  because  they  are  permanently  mounted  to  the 
mirror  in  the  sensor.  In  the  chamber  we  use  for  calibration, 
these  sensors  can  be  as  much  as  0.1°C  away  from  the  cylinder 
temperature.  The  temperature  measured  by  a  sensor  is  found  from 
its  resistance  by  inverting  the  standard  equation 

R  =  Rq(1  +  ot  T)  , 

where  T  is  in  ®C,  Rg  the  resistance  at  the  ice  point,  and  a  is 
the  temperature  coefficient.  For  the  air  temperature  sensors 
a  =  0.003921,*  for  the  dewpoint  sensors  ot  =  0.003895.  The 
calibration  procedure  is  to  determine  the  values  of  RQ  for  each 
sensor  from  the  known  temperature  and  the  measured  resistances. 
This  is  done  at  several  temperatures  to  eliminate  errors  which  can 
be  introduced  by  temperature  inhomogeneities  in  the  chamber. 
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Calibration  accuracy  is  +  0.01°C  for  the  air  and  sea  surface 
temperature  sensors  and  +  0.1  ®C  for  the  dew  point  sensors. 

Wind  Speed;  Calibrations  are  performed  in  one  of  the  NPS  wind 
tunnels  using  a  calibrated  hot  wire  sensor  as  the  standard.  We  do 
not  attempt  to  perform  the  calibration  to  an  accuracy  any  greater 
than  +0.5  m/sec.  We  have  not  found  any  errors  or  differences 
between  sensors  within  0.5  m/sec.  Salt  loading  of  the  bearings  on 
a  sensor  can  cause  it  to  be  outside  of  this  range.  This  is  easily 
detected  long  before  the  calibration  is  affected  by  observing  the 
decay  rate  of  the  shaft  rotation  when  no  cup  is  in  place. 

Affected  sensors  are  never  used. 

Inversion  Height:  The  acoustic  sounder  records  are  calibrated 
with  radiosondes  every  12  hours  during  most  shipboard  operations. 
It  is  misleading,  however,  to  consider  this  to  be  a  good 
calibration  of  the  sounder,  it  is  only  semi-quantitative.  The 
acoustic  return  which  the  sounder  senses  comes  from  reflection 
from  small  scale  temperature  inhomogeneities.  The  inhomogeneities 
are  caused  by  a  combination  of  the  temperature  gradient  and 
mechanical  turbulence  at  the  inversion.  The  return,  thus,  is  not 
necessarily  from  the  height  where  the  temperature  break  occurs, 
which  is  the  inversion  height  deterined  from  the  radiosonde. 
Normally,  radiosonde  and  sounder  determinations  of  boundary  layer 
depth  agree  to  within  40m. 
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